Star Formation at High Angular Resolution 
ASP Conference Series, Vol. S-221, 2003 
M.G. Burton, R. Jayawardhana & T.L. Bourke 



Models of Accretion Disks around Young Stars 



Paola D'Alessio 

Centro de Radioastronomia y Astrofisica, Universidad Nacional 
Autonoma de Mexico, Ap. Postal 3-72, Morelia, Michoacan, Mexico. 



m 
o. 

<^> , Nuria Calvet, Lee Hartmann 



Smithsonian Astrophysical Observatory, 60 Garden Street, MS 42, 
Cambridge, MA 02138 



Oh 

<D 

James Muzerolle 

■ Steward Observatory, University of Arizona, 933 N. Cherry Ave., 
Tucson, AZ 85721 

Michael Sitko 

O . 

Q\ \ Department of Physics, University of Cincinnati, Cincinnati, OH 

m ■ 45221-0011 

o 

^r) [ Abstract. We discuss the importance of accretion in calculating disk 

■ models for young stellar objects. In particular, we show that a disk 
inner rim, irradiated by both the star and the accretion shocks at the 

O stellar surface, can naturally explain recent observations of DG Tau with 

A | the Keck interferometer. We present models for two objects, with mass 

^ ■ accretion rates differing by almost two orders of magnitude, to illustrate 

\ the effects of accretion on the overall disk structure and emission. 

> 

. 1. Introduction 

Disks around Classical T Tauri Stars (CTTS) are accretion disks. This means 
that a fraction of disk material loses energy and angular momentum, falling 
into the gravitational potential well of the star. The standard view of the final 
accretion process in a low mass star is that a well organized stellar magnetic 
field truncates and couples to the inner disk at a few stellar radii (Konigl 1991; 
Shu et al. 1994). Disk material falls onto the star along the stellar magnetic 
field lines in a magnetospheric flow, releasing most of the accretion luminosity, 
L s hock ~ GM t M/fl»(l — R*/Ri) ~ L acc , at the stellar surface in hot accretion 
shocks(Calvet & Gullbring 1998, CG98). The emission of these shocks produce 
a blue/UV continuum excess that veils the stellar photospheric absorption lines. 
Detailed measurement of veiling excesses, combined with models of the hot con- 
tinuum, have resulted in a long list of mass accretion rate estimates in CTTS (see 
Calvet, Hartmann & Strom 2000, and references therein). In addition, models of 
the gas infalling through the magnetosphere have successfully reproduced blue 
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shifted emission peaks, blue line asymmetries and red shifted absorption compo- 
nents of emission lines in CTTS (Hartmann, Hewett & Calvet 1994; Muzerolle, 
Hartmann & Calvet 1998; Muzerolle, Calvet & Hartmann 1998, 2001), providing 
strong support to the magnetospheric model of accretion. 

The average mass accretion rate for CTT disks in several known associations 
is ~ 10~ 8 Moyr" 1 (Gullbring et al. 1998; Hartmann et al 1998.) For this low 
mass accretion rate, the infrared SED produced by the dust in the accretion 
disk is dominated by reprocessed stellar radiation intercepted by the disk curved 
surface (Kenyon & Hartmann 1987; Calvet et al. 2000). Thus, why should one 
get into the trouble of including accretion in a model describing a disk? We 
discuss in this contribution why accretion plays a key role in disk structure and 
cannot be ignored. 

Models of passive disks (i.e., disks heated only by stellar irradiation) have 
been proposed to explain the observed SEDs of CTTS and Herbig Ae stars 
from IR to mm wavelengths (Chiang &: Goldreich 1997). In contrast to power 
law models, where mass surface density and temperature of the disk are given 
by arbitrary power laws, the passive disk models treat in detail the energy bal- 
ance between heating by absorbed stellar radiation and radiative cooling at each 
point. In recent years, powerful radiative transfer techniques have been devel- 
oped for passive disk models, which result in successful accounting for several 
observed features (see Dullemond, van Zadelhoff & Natta 2002; Dullemond 2002; 
Dullemond & Natta 2003). However, even in the most sophisticated versions of 
these models, the disk mass surface density is still assumed to be an arbitrary 
power law of distance, £ = T,q(R/ Rq)~ p , and So and p are adjusted to fit a 
particular observation, with no link to the underlying physics of the disk. 

On the other hand, in a steady accretion disk the mass accretion rate is 
related to the angular momentum flux produced by a viscous torque. As a 
consequence, the mass surface density is given by £ oc M/v t , where v t is the 
viscosity coefficient. It is common to adopt the a prescription (Shakura & Sun- 
yaev 1974), such that v t = ac s H, where c s is the local sound speed, H is the 
local gas scale height and a summarizes the unknown origin of the viscosity. 
Thus, instead of assuming an arbitrary power law for E(i?), the surface density 
is a result of the calculations in accretion disk models, given M (which can be 
independently quantified from observations), the temperature (which is calcu- 
lated from viscous and irradiation heating) and a (which, in principle, can be 
estimated from magnetohydrodynamic simulations; however, in general, is taken 
to be a free constant parameter). 

Irradiated accretion disk models show that even for M ~ 10~ 8 M Q yr _1 , 
viscous dissipation is an important heating mechanism close to the mid-plane 
and close to the central star (R < 1 AU) (D'Alessio et al. 1997, 1999, 2001). 
On the other hand, stellar radiation intercepted by the disk upper layers is the 
most important heating mechanism of the disk atmosphere and the outer disk 
(Kenyon & Hartmann 1987; Calvet et al. 1991, 1992 =CMPD; D'Alessio et 
al. 1998, 1999, 2001). The upper layers of the atmosphere are hotter than 
the interior (CMPD), and this is a characteristic also present in the passive 
models mentioned above. These upper optically thin hot layers are probably 
responsible for some of the observed emission features (silicate 10 [ira band, 
Natta, Meyer & Beckwith 2000; CO overtone emission, Najita et al. 1996). 
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Their reprocessed radiation heats the disk interior, but do not contribute much 
to the continuum, specially when scattering of stellar radiation by atmospheric 
dust grains is included (D'Alessio et al. 1998; Dullemond & Natta 2003). 

In this paper we illustrate the importance of taken into account the accretion 
nature of the disk. First of all, the disk mass surface density is related to its 
mass accretion rate and viscosity, and the local viscous dissipation should be 
included as a heating source, even if the mass accretion rate is small. Second, 
the disk receives radiation from both the star and the accretion shocks at the 
stellar surface. This becomes particularly important for those objects where the 
accretion luminosity is comparable to the stellar luminosity, and has profound 
effects on the structure of their inner disk, as discussed below. It is specially 
relevant for high resolution interferometric observations, and opens up a number 
of possibilities to unveil the physics of accretion. 

2. Irradiated Rim 

Muzerolle et al. (2003a, MCHD) have presented near-infrared spectra of the 
excess continuum emission from the innermost regions of 9 CTTS. The shape 
of this excess is consistent with that of a single temperature black-body with 
T ~ 1200 — 1400 K, which are the expected sublimation temperatures for sili- 
cates under typical disk conditions (Pollack et al. 1994). The observed SEDs 
agree with that of a rim at the dust sublimation radius, receiving radiation 
frontally. A rim irradiated by the central star has been proposed by Natta et 
al. (2001), Dullemond, Dominik & Natta (2001) to explain the observed excess 
at 3 /im in the SEDs of Herbig Ae stars, and by Tuthill, Monnier & Danchi 
(2001) and Monnier Sz Millan-Gabet (2002) to explain interferometric near-IR 
images of disks around Herbig Ae and CTT stars. However, in the CTTS case, 
MCHD find no direct relationship between the radius of the rim and the stel- 
lar luminosity; rather, this radius increases with accretion luminosity. This is 
naturally explained when the fact that the disk is irradiated by both the stellar 
and the accretion luminosity is taking into account (MCHD.) In the following 
section we outline the model for the temperature at the rim atmosphere, which 
can explain the MCHD observations, and we show how the rim radius depends 
on the accretion and stellar luminosities as well as on the dust properties. 

2.1. Model for the irradiated Rim of an accretion disk 

We follow a simplified version of CMPD's treatment to describe the radial tem- 
perature distribution in the rim atmosphere, assuming it has plane-parallel ge- 
ometry and constant opacity. The radiation field is separated into two frequency 
ranges, one characteristic of the disk (and the related quantities have subindex 
"d" ) and one characteristic of the incident radiation (with subindex "i" ) . In the 
case of the rim, the incident field is given by the superposition of the stellar 
radiation field (characterized by the star effective temperature, T*) and the ac- 
cretion shock radiation field (characterized by a temperature ~ 8000 K and a 
luminosity L shock , CG98). 

For simplicity, we assume that the scattering of incident radiation is es- 
sentially forward, which is equivalent to take the albedo to be equal zero in 
CMPD equations (i.e., J s = H s = 0). Also, we assume that the incident radia- 
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tion reaches the rim with an angle relative to the normal to the surface, i.e., 
cos#o = /j,q = 1, and that there is no other heating source than the incident 
radiation. Thus, at every depth into the rim atmosphere the net outward radial 
flux at the disk frequency range is equal to the incident stellar flux, i.e, 

H d (r s ) = ge- (1) 

where Fq = (L* + L s ^ oc fc)/4-7ri? 2 and t s is the optical depth to the stellar radi- 
ation in the incident direction (i.e., along rays parallel to the disk mid-plane, 
perpendicular to a cylindrical surface at the dust sublimation radius). 

The zeroth moment of the transfer equation can be written as (Mihalas 
1978), 

d -^ = J d -S d (2) 
dr d 

where T d is the optical depth at the disk frequency range, S d is the local source 
function and J d is the mean intensity. CMPD assume strict LTE, i.e., S d = 
B d = crT^/ir. Here, we will include a term corresponding to the emissivity 
by scattering at the disk frequency range, i.e., S d = (n d B d + cr d J d )/(n d + cr d ), 
assuming isotropic scattering at this frequency range (Mihalas 1978). With this 
source function, the zeroth moment of the transfer equation is given by 

R _ T XddH d 

D d — J d — • (3J 

K d dT d 

Using the Eddington approximation, the first moment of the transfer equa- 
tion is 

^ = -3H d = (4) 

dT d 47T 

where we have defined q = T s /r d . Integrating J d from eq(4), using J d (0) = 
2H d (0) as the boundary condition, and substituting J d and dH d /dr d from eq(l) 
in eq(3), the temperature as a function of T d can be written as 



T d (r d ) 



qj \ K d q 



(5) 



The temperature at the innermost radius of the rim, where T d = t s = 0, 
should be equal to the silicates sublimation temperature (assuming the inner 
gaseous disk has a negligible radial optical depth). Thus, the rim radius R su b is 



R 



sub 



L* + L s 

hock \ I n Xd 

2 + q — 



167TCTR /V K d 



1/2 i 

jS- ( 6 ) 

sub 



Eq(6) is essentially similar to the equation for the rim radius proposed by 
Tuthill et al. (2001), except for the enhanced heating due to the accretion 
shocks at the stellar surface (given by -L s ^ cfc)) an d the addition of 2 to the 
opacity term inside the parenthesis, which arises from the explicit inclusion of 
the disk radiation field. 

The temperature at the level where r d = 2/3, which characterizes the emer- 
gent flux, is 
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Modeling of the observed near IR SEDs of CTTS indicates that the contri- 
bution of the optically thin layers of the rim (with dust) to the SED is negligi- 
ble, so that the emergent intensity l v differs from B u (T e ff) by only few percent 
(MCDH). Also, we find that T eff is lower than T(0) = T sub within ~ 100 K, so 
that the rim atmosphere is almost radially isothermal (MCDH). 

Finally, the emergent flux is F u = I u fl(i), where is the solid angle 

subtended by the rim as seen by the observer and i is the inclination angle 
between the disk symmetry axis and the line of sight. The rim height is z r i m , 
defined as the height where the optical depth to the stellar radiation impinging 
radially becomes unity at the rim radius R su b- For an accretion disk this height 
can be estimated considering the rim to be vertically isothermal with a mass 
surface density given by £ = M/6irv t [l- {R*/R su b) 1/2 ] (Muzerolle et al. 2003b). 
We find typically 

Zrim ~ 4i^ r j m , where H r i m is the gas vertical scale height of 
the rim (assumed to be vertically isothermal). 

Defining 5 = z r i m tan i / R su b, the solid angle of the rim is 



= it (^j^j 2 cos i for <5>1 (8) 

n(i) =4(^f^j 2 ^( S m- 1 5 + sVT^P) for 5 < 1 (9) 

Both expressions multiplied by B u (T e ff), are consistent with the formulae for 
the emergent flux given in the Appendix B of Dullemond et al. (2001), except 
for our different definition of R su b, which assumes that the dust is optically thin 
at the sublimation radius, and for using a T e f f different from the sublimation 
temperature. 



3. Results 

In this section, we discuss irradiated accretion disk models, including the rim, 
that explain high resolution IR observations or SEDs of different objects. The 
rim model just described fits the near IR spectra of several CTTS, with rim 
radii in the range 0.07-0.54 AU, increasing with higher stellar and accretion 
luminosities (MCDH). Here, we first show that this model can explain the recent 
interferometric observations presented by the Keck team. Then, we present 
models for two systems with widely different accretion luminosities, to illustrate 
the effects of accretion. 

3.1. DG Tau 

Keck interferometric observations of DG Tau (Colavita et al. 2003) showed 
that the system was resolved at K, and derived a size of ~ 0.12-0.24 AU. This 
radius is ~ 3 times larger than the dust sublimation radius inferred assuming 
the rim is irradiated only by the central star. Specifically, for an optically thick 
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wall, R' sub = yL* /16-ko-rT^J . Using typical parameters for a K7-M0 central 

star (T* = 4000 K, R* = 2 R and M* = 0.3 M ), and a silicate sublimation 
temperature T sub = 1400 K, R' sub ~ 0.07 AU. However, DG Tau is a "continuum 
star", with a mass accretion rate M = 5 x 10~ 7 M yr~ x (Gullbring et al. 2000). 
Assuming the typical dust properties derived for the sample of CTTS studied 
by MCHD (i.e., q = 1, Xd/^d = 2), we find that R sub ~ 0.2 AU, which is ~ 3 
time larger than R' sub , and consistent with the observed inner radius of DG Tau. 



3.2. LkCa 15 

We have calculated a model for the structure and emission of LkCa 15, which is 
a low mass accretion rate CTTS. The details of the calculation of an irradiated 
accretion disk model are described by D'Alessio et al. (1998, 1999, 2001). The 
input parameters constrained by observations are M = 1.3 x 10~ 9 M yr -1 , 
T* = 4350 K, R* = 1.5 R , M* = 0.8 M (Hartmann et al. 1998) Figure 1- 
right panel shows the comparison between the observed and predicted SEDs. We 
explain the shape of the observed SED between 60 fim to 1.3 mm, assuming that 
the outer disk (R > 3 AU) has a uniform grain size distribution characterized 
by n(a) oc a~ 3 ' 5 , with a rnax = 1 mm. These are typical dust properties such 
that an irradiated disk model with the median M and a typical central star, 
fits the median SED of CTTS in Taurus ( D'Alessio et al. 2001). We adopt 
a maximum disk radius Rj = 200 AU and a = 0.001, corresponding to a disk 
mass ~ 0.04 M . Qi et al. (2003) find that the mm SED does not show a 
single spectral index, which is naturally explained, by taking into account the 
contribution of optically thick regions at small radii. 

In this model, the disk rim is placed at R su b ~ 0.06 AU. Being > L s h oc k, 
the rim is mostly heated by stellar radiation. There is a shadowed region between 
the rim and ~ 0.16 AU, which produces a negligible contribution to the SED. We 
assume that this region is only being heated by viscous dissipation. Finally, there 
is an intermediate region between ~ 0.16 and 3 AU, where the disk atmosphere is 
characterized by smaller grains than the disk interior and the disk outer regions. 
This region is responsible for the observed 10 \xm silicate band which is formed 
in the upper layer, heated by direct stellar radiation (Natta et al. 2000). In 
order to fit the shape of the silicate band, we have used optical constants for 
olivine, MgFe Si04 (Jager et al. 1998). To produce a continuum around 10 
fim consistent with the observed value, the dust to gas mass ratio of the upper 
atmosphere should be a factor ~ 0.05 smaller that the standard value, probably 
reflecting that dust growth and settling has taken place in the inner disk. In 
principle, the details of this model can be tested by high resolution observations 
capable of resolving scales between 0.4 and 20 mas. 



3.3. DR Tau 

DR Tau is a continuum star, characterized by a strong UV/optical continuum. 
We have calculated a model for the disk of DR Tau, with the following parame- 
ters (MCHD) T* = 4000 K, i?* = 1.9 R Q and M* = 0.9 M for the central star. 
We take M = 5 x 10" 7 M yr 1 , roughly consistent with the value inferred from 
the optical/UV excess (Kenyon et al. 1994; Gullbring et al. 2000). A uniform 
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Figure 1. Right panel: SED of LkCal5 showing different contribu- 
tions. Left panel: SED of DR Tau showing the contributions of differ- 
ent zones. For comparison, also a SED of the same disk model only 
irradiated by the central star is shown. See D'Alessio et al. (2003) for 
references to observational data. 

maximum grain size a max = 1 mm, a = 0.01, = 50 AU, corresponding to a 
disk mass ~ 0.2 M . In contrast to the case of LkCa 15, accretion luminos- 
ity determines the structure and emission of the entire disk for this object. The 
rim radius is R su b ~ 0.15 AU, and it is heated mostly by the hot continuum 
produced in the accretion shocks at the stellar surface. Figure 1-left panel shows 
that the rim emission dominates the SED between 3 and 6 fim. Local viscous 
dissipation is important in heating the disk photosphere for R <J 3 AU and heat- 
ing the mid-plane at every radius (R <J 50 AU). However, the photosphere of 
the outer disk is heated mostly by radiation from the accretion shocks. This 
model can also be tested by high resolution observations. 

4. Conclusions 

We have shown that accretion plays an important role in determining the struc- 
ture and emission of disks around CTTS and cannot be dismissed. The surface 
density of the disk depends on the mass accretion rate, and moreover, accretion 
energy determines the location of the transition between dust and gas in the 
inner disks, with direct implications for interferometric observations. In partic- 
ular, by including the contribution of accretion, the puzzling Keck interferometer 
observations of DG Tau, with an disk inner radius larger than the sublimation 
radius expected from stellar luminosity alone, can be understood. 
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